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Progress  in  a  number  of  ongoing  projects  is  reported.  These  include 
studies  of  bistatic  and  monostatic  surface  scattering  strength,  surface 
reverberation  modelling,  measurement  of  slope  and  correlation  statistics  of 
surfaces,  and  measurement  of  Doppler  shifts  and  spread. 

The  extension  of  previous  scattering  strength  measurements  to  lower 
frequencies  has  shown  the  strengths  observed  to  be  10  to  15  dB  higher  under 
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similar  wind  conditions  but  with  wavelengths  approximately  8  times  longer. 
These  studies  have  led  to  the  beginning  of  development  of  a  model  for 
reverberation  based  on  a  characterization  of  the  surface  as  an  ensemble  of 
facets. 


A  method  for  precision  measurement  of  model  tank  waveheight  and  wave 
slope  time  series  has  been  developed.  The  system  is  based  on  fine  wire 
capacitance  probes.  It  has  been  used  to  compile  a  set  of  slope  measure¬ 
ments  which  is  consistent  with  data  reported  by  other  investigators.  A 
system  for  automatic  measurement  of  time-space  waveheight  correlation  has 
been  developed. 

A  seventeen  element,  320  kHz.  linear  array  having  a  beam  width  of  4° 
and  sidelobe  ejection  of  23  dB.  has  been  designed  and  constructed.  The 
array  uses  foam  backed  PZT  elements  and  is  intended  for  low  grazing  angle 
measurements . 

Additional  measurements  of  the  alignment  and  wind  direction  depen¬ 
dence  of  frequency  spread  and  mean  shift  in  forward  scatter  have  been 
made.  These  measurements  indicate  that  shift  and  spread  are  quite  sensi¬ 
tive  to  alignment  error  in  crosswind  orientations  but  that  even  under 
conditions  of  precise  alignment,  some  frequency  shift  does  occur. 
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!•  In troduc t ion :  During  the  contract  year  as  amended  a  number  of 
projects  were  initiated  or  continued  in  support  of  the  general  con¬ 
tract  objective,  characterization  of  acoustic  scattering  strength  from 
wind  roughened  water  surfaces.  This  dociment  is  a  statement  of  pro¬ 
gress  as  of  the  end  of  the  contract  period  for  each  of  these  projects. 

The  principal  effort  of  the  program  has  been  the  continuation  of  an 
empirical  study  of  the  scattering  strength  of  wind  driven  surfaces 
under  bistatic  as  well  as  monostatic  configurations.  This  study  was 
initiated  to  provide  a  characterization  of  this  statistic  in  the 
absence  of  any  Marking  theoretical  formulation  so  as  to  permit  the 
refinement  of  global  propagation  models  to  include  surface  reverbera¬ 
tion  effects.  During  previous  contract  years  a  precision  goniometer 
was  constructed  to  permit  efficient  measurement  of  scattering  strength 
under  a  variety  of  geometries  and  an  extensive  study  of  scattering 
strength  at  short  wavelengths  was  conducted  iZornig,  19781.  During  the 
past  year  the  goniometer  has  been  improved  and  a  series  of  measure- 


roents  ot  scattering  strength  at  much  longer  wavelengths  has  been 
condts'ted.  These  measurements  have  revealed  some  substantial  dlifer- 
enees  In  the  behavior  ot  t  tie  statistic  at  these  wavelengths  which  will 
be  described  later  In  tills  report  and  which  Indicate  the  need  for 
t  ur  t  he  r  s  t  ud  v  . 

Motivated  In  part  by  the  results  of  the  empirical  studies  at  this 

facility  and  In  part  hv  the  work  ot  other  invest  iga  tot  s ,  a  basis  tot  a 

theoretical  model  has  been  tormulated  during  the  past  tew  months  which 
shows  promise  of  explaining  some  or  all  ot  the  results  ot  the  experi¬ 
mental  program.  Hits  wvrk,  which  is  an  extension  ot  previous  work  ot 

Patterson  list's]  and  Novarlnt  and  Medwin  attempts  to  explain 

the  reverberant  field  ot  an  lnsonitled  rough  surface  as  the  result 
dtttractlve  Interactions  with  an  Idoallrcd  piecewise  plane  surface  the 
statistics  ot  which  can  be  derived  trom  measurements  ot  the  height  and 
•lope  statistics  of  the  actual  surface.  The  model,  while  bv  no  means 
fully  developed,  does  Incorporate  effects  which  are  cleat Iv  important 
In  reverberation  at  low  grarlng  angles  and  long  wavelengths  but  which 
are  neglected  In  Helmholtr  Integral  based  t ormul at  tons . 

In  a  separate  studv  also  continuing  previously  reported  work,  a  numbet 
of  additional  measurements  of  floppier  spread  In  the  high  roughness 
domain  have  been  made.  These  measurements  wv  i  e  made  In  order  to  clar¬ 
ify  the  relationship  between  beam  pattern  orientation,  orientation 
with  respect  to  wind  direction,  and  mean  floppier  shift.  Hie  results  ot 
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this  series  of  experiments,  reported  later  In  this  docunent,  Indicate 
that  the  prediction  of  Doppler  shift  and  spread  will  probably  require 
a  more  precise  model  than  those  presently  available.  This  work  has 
motivated  a  theoretical  Investigation  supported  under  other  auspices 
which  will  be  published  in  the  near  future  [Tuteur  et  al.,  1979]. 

The  acoustic  studies  conducted  during  previous  contract  years  have  all 
indicated  that  a  more  thorough  characterization  of  the  scattering 
surface  would  be  required  If  any  general  model  capable  of  explaining 
the  phenomena  observed  were  to  be  developed.  Accordingly  studies  of 
the  surface  slope  statistics  and  time-space  covariance  function  were 
undertaken  and  will  be  documented  In  this  report.  These  measurements 
were  made  using  a  refinement  of  a  uavehelght  measuring  apparatus  which 
has  been  developed  over  the  past  several  years  based  on  a  capacitance 
wire  detection  scheme.  They  exemplify  a  very  high  precision  vaveheight 
measurement  capability  which  will  used  to  extend  the  documentation  of 
water  surfaces  In  the  future. 

During  past  discussions  with  project  sponsors  of  the  direction  of  the 
modelling  project  at  this  facility,  a  recurrent  theme  has  been  the 
issue  of  very  low  grazing  angle  measurement  of  scattering  statistics. 
During  the  contract  period  covered  by  this  report  funds  were  requested 
for  partial  support  of  the  development  of  an  ultrasonic  array  capable 
of  producing  a  narrow  sidelobeless  beam  pattern  for  use  in  future 
studies  of  heavily  shadowed  scattering  processes.  The  array  was  de- 
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signed,  constructed  and  characterized  during  the  past  year  with  suc¬ 
cess  and  Is  described  In  this  report. 

2.  B1 s t a t lc  Scattering  St rength  St ud les :  This  section  describes  the 
results  of  a  series  of  measurements  of  the  scattering  strength  of 
short  fetch  model  tank  surfaces  under  a  variety  of  experimental  situ¬ 
ations.  These  measurements  were  conducted  using  a  technique  described 
previously  (Zornig,  1978]  for  the  case  of  short  wavelengths.  The 
acoustic  wavelengths  used  in  this  case,  however,  were  much  longer  than 
those  of  the  earlier  study.  In  addition,  scattering  strengths  were  not 
measured  for  bistatic  differential  azimuths  greater  than  105  deg. 

(.near  forward).  These  measurements  were  omitted  because  of  difficulty 
in  obtaining  data  in  near  forward  orientations  that  were  free  of 
direct  path  contributions  and  because  the  conceptual  basis  of  the 
statistic  is  of  doubtful  usefulness  in  near-forward  orientations  in 
which  contributions  to  the  scattered  field  cannot  be  assumed  to  be 
uniformly  distributed  over  the  area  of  the  surface  insonified.  These 
orientations  are  difficult  to  characterize  primarily  because  of  the 
absence  of  a  simple  statistic  that  is  Independent  of  the  projector/ 
hydrophone  characteristics. 


The  procedure  used  in  these  experiments  remains  essentially  unchanged. 
A  directional  projector  is  used  to  insonify  a  portion  of  the  wind 
driven  surface  from  a  specified  direction  at  a  specified  wavelength  at 


a  relatively  well  defined  frequency  (more  than  20  cycles  per  pulse). 
The  projector  used  was  a  piezoelectric  transducer  having  an  effective 
aperture  of  7.7  cm.  and  a  principal  resonance  of  90  kHz.  The  transdu¬ 
cer  exhibited  a  reasonably  piston-like  beam  pattern  at  the  principal 
resonance  and  had  sufficient  response  to  be  useful  at  180  kHz  and  36u 
kHz  as  taell.  The  beam  patterns  at  these  harmonic  frequencies  appeared 
to  be  piston-like  as  well,  but  were  somewhat  less  ideal  than  that 
measured  at  90  kHz.  A  directional  hydrophone  is  used  to  sense  the 
scattered  field  at  some  other  orientation  during  a  precisely  defined 
and  constant  time  Interval.  The  combination  of  beam  patterns  and  time 
gate  are  used  to  compute,  for  each  orientation,  an  effective  insoni- 
fication  pattern  on  the  surface.  This  pattern  together  with  a  cali¬ 
bration  measurement  made  w>ith  the  same  apparatus  from  a  smooth  surface 
in  a  specular  reflection  orientation  is  used  to  calculate  from  the  raw 
measured  returns,  the  scattered  power  that  would  have  been  received  at 
one  yard  from  a  uniformly  insonlfied  square  yard  of  surface  due  to 
unit  incident  power.  This  statistic  is  called  the  scattering  strength 
of  the  surface.  The  calculation  of  the  statistic  also  depends  on  the 
spreading  law  governing  reradiation  from  the  surface.  However,  the 
experimental  geometry  was  such  that  an  error  in  the  assumed  form  of 
the  spreading  law  would  have  made  only  a  small  difference  in  computed 
scattering  strength.  There  has  been  and  still  remains  considerable 
controversy  over  the  law  governing  spreading  from  rough  scattering 
surfaces  [Mikeska  and  McKinney,  1978]. 
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Using  this  procedure  scattering  strengths  were  measured  for  a  nimber 
of  differential  azimuths  (the  difference  between  projector  and  recei¬ 
ver  orientations  with  respect  to  the  center  of  the  scattering  area)  at 
each  of  tv*>  wind  speeds,  tw>  grazing  angles,  and  three  projector 
orientations  with  respect  to  the  wind  (upwind,  crosswind,  downwind). 
The  surface  statistics  are  the  same  as  those  used  previously  and  are 
further  documented  later  in  this  report.  This  basic  set  of  data  is 
presented  as  Figures  1  -  12.  The  plots  differ  from  those  presented  in 
previous  reports  in  that  the  backscatter  direction  is  at  the  left  end 
of  the  abscissa.  From  these  data  two  basic  observations  are  apparent: 

a.  Scattering  strengths  are  apparently  substantially  higher  at 
this  longer  wavelength.  Figure  13.  shows  the  scattering  strength 
for  backscatter  from  the  Wind  4  surface  at  three  projector 
orientations,  two  grazing  angles  and  four  wavelengths.  It  is 
clear  from  these  results  that  the  change  of  wavelength  has  pro¬ 
duced  what  appears  to  be  a  consistent  10  to  15  dB.  increase  in 
scattering  strength  at  all  orientations.  The  increase  is  appar¬ 
ently  higher  at  higher  grazing  angles,  but  whether  it  is  propor¬ 
tional  or  not  is  an  open  question.  As  will  be  seen  below,  a 
recurring  theme  in  this  report  is  the  speculation  that  reverber¬ 
ation  at  low  grazing  angles  is  primarily  due  to  diffractive 
effects.  If  we  believe  that  the  surface  behaves  as  an  ensemble  of 
plane  facets,  each  reradiatlng  with  a  pattern  that  is  dependent 
on  its  size  and  the  incident  wavelength,  then  it  seems  reasonable 
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Figure  3-  -  Scattering  strength  measured  at 
177  kHz.,  12  deg.  grazing  angle,  projector 
facing  upwind,  wind  speed  8.3  m./sec. 
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Figure  4.  -  Scattering  strength  measured  at 
177  kHz.,  1 2  deg.  grazing  angle,  projector 
facing  downwind,  wind  speed  5.4  m./sec. 
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Figure  5.  -  Scattering  strength  measured  at 
177  kHz.,  12  deg.  grazing  angle,  projector 
facing  crosswind,  wind  speed  m./sec. 


Figure  6.  -  Scattering  strength  measured  at 
177  kHz.,  12  deg.  grazing  angle,  projector 
facing  upwind,  wind  speed  5.^  m./sec. 
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Figure  7.  -  Scattering  strength  measured  at 
177  kHz.,  30  deg.  grazing  angle,  projector 
facing  downwind,  wind  speed  8.3  m./sec. 
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Figure  8.  -  Scattering  strength  measured  a 
177  kHz.,  30  deg.  grazing  angle,  projector 
facing  crosswind,  wind  speed  8.3  n./sec. 
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Figure  9.  -  Scattering  strength  measured  at 
177  kHz.,  30  deg.  grazing  angle,  projector 
facing  upwind,  wind  speed  8.3  m./sec. 
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that  at  longer  wavelengths  the  patterns  would  become  more  omni¬ 
directional  and  would  tend  to  cast  more  power  in  non-specular 
directions.  However,  the  degree  to  which  this  is  ttie  case  wwuld 
depend  on  the  sizes  oi  such  facets  and  the  distribution  oi  their 
occurance . 

b.  At  this  wavelength  the  dependence  of  scattering  strength  on 
projector  orientation  with  respect  to  the  wind  remains  the  same. 
It  was  noted  during  the  short  wavelength  series  of  experiments 
that  a  difference  of  as  much  as  15  dB.  in  scattering  strength 
could  be  observed  between  upwind  and  downwind  backscattering 
strengths.  This  was  supposed  to  be  due  to  the  asymmetry  of  the 
short  fetch  wind-driven  surface  used  in  these  experiments.  One 
could  speculate  that  at  some  wavelength  the  effect  of  this  asym¬ 
metry  might  diminish,  but  at  this  point  there  is  apparently  no 
change . 

In  addition  to  the  basic  series  described  above,  several  additional 
scans  were  made  at  wavelengths  twice  and  one  half  as  large.  These  runs 
were  exploratory  in  the  sense  that  no  extensive  scan  at  these  wave¬ 
lengths  was  contemplated  unless  anomalous  behavior  was  observed.  They 
are  given  as  Figures  14  -  19.  As  can  be  seen  from  a  comparison  of 
Figures  14  and  15  with  Figures  10  and  11,  reducing  the  wavelength  by  a 
factor  of  two  had  only  the  effect  of  reducing  the  scattering  strength 


at  both  crosswind  and  downwind  orientations.  On  the  other  hand  doub- 
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Figure  15.  -  Scattering  strength  measured  at 
92  kHz.,  30  deg.  grazing  angle,  projector 
facing  crosswind,  wind  speed  8.3  m./sec. 
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Figure  16.  -  Scattering  strength  measured  at 
357  kHz.,  30  deg.  grazing  angle,  projector 
facing  downwind,  wind  speed  8.3  m./sec. 
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Figure  18.  -  Scattering  strength  measured  at 
302  kHz.,  30  deg.  grazing  angle,  projector 
facing  downwind,  wind  speed  8.3  m./sec. 
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Figure  19.  -  Scattering  strength  measured  ; 
302  kHz.,  30  deg.  grazing  angle,  projector 
facing  crosswind,  wind  speed  8.3  m./sec. 
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ling  the  wavelength  (Fig.  16,17)  seems  to  have  had  the  effect  of  both 
raising  the  scattering  strength  and  reducing  the  downwind/crosswind 
difference.  Thl6  observation  is  admittedly  based  on  an  anomaly  in  just 
one  run,  but  ft  does  suggest  that  further  exploration  at  even  longer 
wavelengths  may  be  worthwhile. 


3.  The  Facet  Model :  During  the  past  year,  in  response  to  the  continu¬ 
ing  need  for  a  viable  theoretical  framework  on  which  to  hang  the 
considerable  mass  of  data  that  has  been  collected,  a  study  has  been 
conducted  to  determine  the  possibility  of  constructing  a  semiempir ical 
model  for  reverberant  scattering  based  on  the  notion  of  the  surface  as 
an  ensemble  of  facet-like  scatterers.  The  notion  of  modelling  a  scat¬ 
tering  surface  as  an  ensemble  of  statistically  defined  scatterers  has 
been  championed  for  some  years  by  Middleton  [1977],  but  has  not  been 
used  widely  due  in  some  part  to  the  difficulty  of  characterizing  the 
scatterers  in  terms  of  measurable  characteristics  of  the  physical 
surface.  Patterson  [1963,1964]  and,  more  recently,  Novarini  and  Medwin 
[1978]  have  presented  limited  models  of  the  surface  based  on  the 
notion  of  plane  facets,  but  these  models  are  relatively  primitive  and 
do  not  in  their  present  condition  permit  the  quantitative  prediction 
of  experimental  data. 


The  basic  approach  of  the  facet  model  is  to  represent  the  surface  as  a 
connected  set  of  plane  reflectors  of  a  specific  size  and  shape.  In 
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that  way  the  surface  would  be  simply  parameterized  and  much  of  its 
complexity  would  be  eliminated.  The  connection  to  the  pnyslcal  surface 


is  made  by  measurement  of  surface  slope  and  other  surface  statistics, 
permitting  the  selection  of  a  distribution  of  facet  sizes  and  slopes 
in  the  model.  The  scattered  field  is  then  expressed  as  the  coherent 
sum  of  the  waves  reradiated  from  each  facet  to  the  measurement  point. 
Each  facet  is  treated  separately  as  regards  its  response  to  incident 
field  and  the  coherent  summation  of  reradiated  fields  is  relied  upon 
to  account  for  the  dependence  between  facets.  The  model  is  clearly 
approximate  and  therefore  lacks  the  elegance  of  the  solution  of  a 
scattering  integral  approach.  On  the  other  hand  formulations  based  on 
such  a  model  can  be  developed  in  stages  of  complexity  as  needed  to 
arrive  at  a  practical  predictor.  The  principal  appeal  of  the  facet 
model  approach  is  that  it  preserves  to  a  large  extent  those  scattering 
phenomena  which  are  assumed  away  in  Helmholtz  integral  formulations 
and  which  show  greatest  promise  in  the  explanation  of  low  grazing 
angle  reverberation  (i.e.  diffraction). 

The  first  step  in  the  development  of  the  facet  model  has  been  the 
derivation  of  an  expression  for  the  contribution  to  the  scattered 
field  of  one  member  of  the  ensemble.  The  simplest  shape  for  that 
single  facet  is  a  rectangle  that  is  very  long  in  one  dimension.  A 
suitable  coordinate  system  is  shown  in  Fig.  20.  S  and  R  refer  to  the 
source  and  receiver  locations  in  the  Z-X  plane  and  a  and  p  are  the 
angles  from  the  Z  axis  to  the  source  and  receiver  directions  respec- 


dimensional  facet  interaction. 


» 
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tively.  Tlie  facet  Is  assumed  to  extend  a  threat  many  wavelengths  In  the 
Y  direction.  The  source  is  assimed  to  be  at  a  distance  sufficient  so 
that  the  insonifying  wavefront  over  the  X  extent  of  the  facet  is 
essentially  plane.  This  is  not  a  severe  restriction  given  that  facets 
themselves  can  be  expected  to  be  relatively  small.  The  amplitude, 
phase,  and  direction  of  the  insonifying  wave  is  determined  by  the 
projector  orientation  and  beam  pattern. 

The  reradiated  field  is  calculated  under  these  assumptions  for  the 
facet  alone.  The  received  pressure  is  found  to  be  [Kerr,  1951;  Good¬ 
man,  19c>8;  Patterson,  1963,41 

p^  cos  6  sin(S  kX) (sin  a  -  sin  t?) 

P(a.p)  “  (1) 

•  exp(lk(K^  +  K,))/S  R,(sin  a  -  sin  t>) 

where  k  is  the  the  acoustic  wavenumber,  X  is  the  facet  size,  Kj  and  R., 
are  the  facet  center  to  source  and  receiver  distances  respectively  and 
p^  is  the  pressure  incident  on  the  surface.  In  the  forward  specular 
direction  a  ■  P  and  the  pressure  PFS  is 

PFS  "  P^  kX  exp(ik(Rj  +  R2>)  cos  a/R,  (2) 

In  the  backscatter  direction  a  »-8and  the  pressure  P^s  is 


% 


P_c  “  P  cos  a  sin(kX  sin  a)exp(ik(R,  +  R,))/sin  a  (3) 

Ob  o  I 


♦ 


1 1 


For  large  values  of  kX  the  received  amplitude  can  be  seen  to  be  very 
small  unless  a  and  j3  are  nearly  equal  (specular  reflection).  As  kX 
becomes  smaller  the  facet  tends  to  scatter  more  omnidirectionally, 
approaching  a  point  scatterer.  In  that  case  the  scattered  intensity 
approaches 


|p  |  ^  cos‘a  (kX)"- 
o 


The  eventual  goal  of  a  facet  model,  of  course,  is  to  predict  the 
scattered  intensity  due  to  a  continuous  rough  surface.  The  model 
represents  that  surface  as  an  approximation  in  terms  of  plane  facets 
parameterized,  in  this  case,  by  size,  slope  and  position.  Addition¬ 
ally,  the  facets  are  constrained  to  be  joined  at  the  edges.  In  the 
case  of  a  random  rough  surface  the  parameters  are  random  variables, 
not  necessarily  independent.  The  scattered  field  due  to  this  surface 
is  calculated  by  summing  the  received  fields  due  to  each  facet  assum¬ 
ing  that  each  is  insonified  by  a  plane  wave  of  amplitude  a^  and  phase 
determined  by  the  source  beam  pattern.  This  way  of  accounting  for 
the  characteristics  of  the  source  is  similar  to  the  approach  used  bv 
Clay  and  Medwin  [1970]  to  include  the  Fresnel  approximation  into  the 
Helmholtz  iniegral  based  formulation. 


(4) 


j 


In  the  case  of  a  random  surface  the  basis  for  parameter  assignment  in 
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the  model  will  be  measurements  of  surface  spatial  statistics.  From 
these  a  model  surface  can  be  generated  (Fig.  21).  In  the  case  of  a  one 
dimensional  representation,  the  total  received  amplitude  due  to  a 
finite  insonified  region  of  the  surface  having  length  L  can  then  be 
expressed  as  the  sum  of  contributions  from  the  facets. 


P  “  I  exp(ik(rj1+Xj2)) (cos  a^  cos  Bj) 


sin  *skX  •  (sin  a^  •  sin  Bj)/*s(sin  -  sin  Bj) 


The  received  intensity  is  then 


1  i  2  2 

|PP*|  -  — —  \  Z  a  cos  a  cos  B.  sin  [  (kX .  /2)(sin  a,  -  sin  B.)] 
R2Z  f  j-1  3  3  3  3  3 


/(sin  Oj  -  sin  B^)' 


+  8  Re  1  *  \  ajam  exp[ikCr  +  r  -(r  -  r  ))] 

'  j»l  m*l  J  J  J 


(cos  a.  cos  B.  cos  a  sin  6  )  sin[  — t*-  (sin  a.  -  sin  B.)] 
j  J  m  m  2  j  j 


sin  [  — —  (sin  a  -  sin  8  )]  /  (sin  a.  -  sin  6.)  (sin  a  -  sin  6); 

2  m  m  j  j  m  m| 


The  constraints  on  the  generation  of  the  surface  are  expressed  as 
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I  X  cos  6.  ■  L 

j-1 


ZJ  +  Xj  Sl"  6J  ‘  Zj+2  '  Xj+1  stn  °j+l 


4J  +  XJ  C°S  l>j  “  4J+1  '  Xj+1  COS  °j+l 


where  z^  Is  the  vertical  displacement  of  facet  center  i  and  ^  is  the 
horizontal  displacement  of  the  facet.  The  facet  slope  is  tarn',.  The 
first  term  in  Eq  .  b  represents  the  sum  of  received  intensities  from 
each  facet  separately.  The  second  term  contains  the  phase  dependent 
interference  contribution  (constructive  or  destructive)  among  facets. 


Let  us  consider  the  case  of  a  distant  projector  in  which  the  insoni- 
fying  wavefront  is  plane.  The  angles  Oj  and  p,  are  expressed  in  terms 
of  the  source  and  receiver  grazing  angles,  }  and  and  the  facet 
slope  ,6  .  ,  by 


'n  -  'i  -  ej 


it/2  -  Y,  +  6. 


For  purposes  of  presenting  a  simple  example  we  further  stipulate  that 
the  source  and  receiver  are  at  the  same  range,  R,  and  grazing  angle, 
>.  Under  these  conditions,  the  intensity  measured  in  the  specular 
direction  (forward  scatter)  is 


(7) 

(8) 

(9) 


(10) 
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L  -  — .  I  [tan2(y-0  )sin2(kX  cos (>-6  ))]  + 

R,  !  j-1  J  J  J 


i  j-1 
2  1  Z 

j-2  m«l 


cos  4.  sin(y-0 . ) sin(y-0  )sin(kX.  cos(y-0  )sin(kX  cos(y-0  )) 
Jm _ j _ J _ j _ 5 _ E _ 


cos(y-0 .)cos(y-0  ) 
j  m 


(ID 


and  the  monostatic  backscattered  intensity  is 


|2  (  i  sin(y-0.)sin(Y+© ,)sin2 (kX.  sin  y  sin  0.) 

.  i  f - J - J - J - J-  ]  + 

o  l  **  L  0  0  J 


‘“S  R,2  h-1 


sin2  y  sin2  0, 


£  j-1 

2  1  I 

j«2  m=l 


cos  4!  (sin(y-0.)sin(y+6,)sin(y+0  )sin(y-0  )) 
_ 15 _ j  j _ 5 _ m _ 


sin  y  sin  0,  sin  0 

j  m 


sin(kX,  sin  y  sin  0.)  sin(kX  sin  y  sin  0  ) 

j  j  m  m 


(12) 


The  term  4  is  the  phase  difference  due  to  path  length  between  facets 


in  backscatter  and  is  given  by 


m-1 


4jm  “  k{Xj  C08^-6j>  +  Xm  cos(Y-6m>  +  2  n^j+1  Xn  co«(yJ-0n>}  (13) 


Similarly  the  term  4'^  for  the  specular  direction  is  given  by 


N 


The  phase  term  i'  .  becomes  very  small  when  y  or  6  and  k  are  small, 
jm 

This  occurs  in  the  case  of  relatively  low  surface  roughness,  a  vir¬ 
tually  smooth  surface.  In  that  case  the  expressions  reduce  to  a  Huy- 


gen's  source  representation  of  a  smooth  surface. 


In  the  case  of  monostatic  backscatter,  the  phase  term  depends  on 
the  distance  between  facet  centers  along  the  direction  of  propagation 
of  the  incident  wave.  In  the  case  of  small  slopes  and  low  grazing 
angles  this  is  approximately  the  horizontal  distance  between  facets. 

If  the  distance  between  scatterers  is  typically  much  smaller  than  a 
wavelength,  then  one  would  expect  that  the  contributions  from  facets 
would  be  coherent  within  Fresnel  zones.  Also,  if  the  facet  sizes  are 
much  smaller  than  a  wavelength,  then  they  behave  as  more  or  less 
omnidirectional  scatterers.  Under  these  conditions  one  would  expect 
that  the  contributions  of  many  Fresnel  zones  would  be  roughly  of  equal 
amplitude  and  uniformly  distributed  phase  and  would  interfere  at  the 
receiver.  At  the  other  extreme,  short  wavelengths  and  relatively  large 
facets,  one  would  ..xpect  largely  specular  reflections  from  facets 
adding  in  a  way  that  is  dependent  on  geometry  and  surface  statistics, 
but  again  very  little  backscatter.  In  the  transition  region,  however, 
it  seems  likely  that  there  exists  a  regime  in  which  wavelengths  are 
long  enough  to  give  rather  broad  reradiation  patterns,  but  in  which 


facet  contributions  add  somewhat  coherently 


1<> 


This  section  hns  presented  a  preliminary  torn  of  a  model  based  on  a 
repr esentat ion  of  a  rough  surface  as  an  ensemble  of  plane  segments.  In 
order  to  be  effective,  we  expect  that  the  model  will  have  to  be  re¬ 
fined  and  extended  to  2  dimensions.  While  some  of  this  may  be  accom¬ 
plished  analytically,  it  is  expected  that  the  end  product  will  be  a 
partially  numerical  model,  because  of  the  simple  parameterization  of 
the  Initial  form  of  the  model,  we  expect  that  numerical  experimenta¬ 
tion  will  permit  development  of  the  model  to  be  guided  by  a  more  or 
less  cont inuous  experimental  validation  process. 


4.  Wa  ter  Surface  Statistics:  Tills  program  of  measurements  has  been 
motivated  by  the  rather  demanding  requirements  of  current  analytical 
models  of  surface  scatter  for  exact  knowledge  of  the  statistics  of  the 
scattering  surface.  One  of  the  principal  advantages  of  model  tank 
studies  Is  that  they  permit  the  use  of  a  repeatable  calibrated  sur¬ 
face.  This  advantage  Is  greatly  reduced  If  that  calibration  is  only 
partial  or  Is  Imprecise.  However,  measurement  of  the  behavior  of  a 
water  surface  in  which  there  are  significant  features  with  sizes  on 
the  order  of  a  millimeter  is  difficult  and  has  been  an  ongoing  effort 
at  this  facility  since  Its  establishment.  A  number  of  current  models 
of  underwater  acoustic  scattering  and  frequency  spreading  depend  on  a 
knowledge  ot  surface  height  and  surface  slope  distributions  as  well  as 
the  second  order  spatial  and  temporal  correlation  functions.  It  hns 
been  suggested,  for  example,  that  backscatter  Is  due  to  reflection 


f row  portions  ot  tho  surtace  having  very  stoop  slopes,  and  that  tor 
similar  reasons,  observed  asvmmottv  In  llopplei  spread  might  ho  duo  to 
asymmetry  in  tho  surface  slope  probability  densities.  The  so  sugges¬ 
tions,  as  veil  as  tho  requirements  ot  tho  facet  model  described  above, 
have  generated  a  requirement  tor  a  detailed  knowledge  ot  surtace  slope 
statistics  in  the  model  tank  facility.  Such  measurements  have  in  the 
past  been  made  by  opt  teal  means  [llobb  et  al .  ,  1*171;  Cox  and  Munk , 

1118;  Cox,  11S8;  Schoolev,  11(<0;  Sturm  and  Sorrell,  11?  1;  Wu ,  1171, 
1177),  but  data  ot  ihls  sort  lacks  both  the  precision  and  statistical 
sufficiency  tor  out  purposes.  This  section  describes  the  development 
of  an  improved  measurement  Instrumentation  and  some  recent  measure¬ 
ments  ot  water  surtace  slopes  and  time-space  covariance  functions. 

An  alternative  to  optical  measurement  is  the  use  ot  wive  staffs,  lit 
the  case  ot  the  small  scale  measurements  used  in  model  tank  experi¬ 
mentation,  the  staft  is  a  very  fine  wire  and  is  normally  called  a 
probe  wive.  The  probe  wire  penetrates  the  surtace  vertically  and  Is 
used  to  sense  the  height  ot  the  watei  at  that  location  as  a  (unction 
oi  time.  These  time  series  ol  wove  heights  contain  information  not 
derivable  irom  glitter  photographic  measurements  and  are  the  basis  lot 
the  data  contained  in  this  section. 

As  before,  tho  water  surface  measurements  were  made  in  the  8  meter 
diameter  by  1.1  meter  deep  model  tank  at  tlx-  Yale  facility.  The  wind 
driven  waves  which  are  the  sub  led  ot  the  measurements  at  e  generated 
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by  a  collimated  fan  blowing  over  the  water  with  approximately  3  meters 
of  fetch.  The  air  velocity  is  regulated  by  a  combination  of  lan  speed 
and  aperture  adjustment.  The  formation  of  standing  waves  in  the  appa¬ 
ratus  is  suppressed  to  a  considerable  extent  by  a  downwind  beach  and 
by  irregularities  in  the  skirts  of  the  tunnel  used  to  contain  the 
wind.  Wind  speeds  reported  wore  determined  by  an  anemometer  suspended 
about  18  cm.  above  the  water  surface  near  the  measurement  site. 

Water  height  measurements  were  made  by  measuring  the  capacitance 
between  the  conductor  of  a  .07  mm.  uniformly  insulated  (Kormvar  in¬ 
sulated  magnet)  wire  and  the  water  in  which  it  is  suspended  vertical¬ 
ly.  The  capacitance  of  such  a  system  is  linearly  proportional  to  the 
immersion  depth  providing  that  the  wire  and  insulating  coating  are 
uniform.  In  order  to  develop  useful  time  series  of  writer  depth,  t  he 
capacitance  must  be  sampled  at  a  rate  on  the  order  of  several  hundred 
times  per  second.  This  was  accomplished  by  using  the  probe  as  the 
capacitor  in  a  high  pass  tiller  with  the  water  at  ground  potential  as 
shown  in  Figures  22  and  23.  Measurements  of  surface  slopes  were  made 
using  tvxi  probes  with  Independent  electronics.  A  constant  amplitude  S 
kHz.  sinusoid  is  applied  to  the  basic  probe  circuit.  The  amplified 
output  of  the  probe,  an  amplitude  modulated  signal,  is  then  detected 
to  derive  a  waveheight  signal.  A  signal  proportional  to  slope  is 
developed  by  subtracting  the  detector  outputs  of  two  piobes  of  equal 
calibration.  Typically  the  result  is  digitized  at  a  rate  of  100  times 
per  second  to  a  precision  of  12  bits  and  stored  for  analysis. 


Figure  22.  -  The  two  probe  waveheight  measure 
ment  system  For  acquiring  time  series  of 
'Jirectional  slopes. 


Figure  23.  -  Circuit  diagram  of  the 
detector.  Calibration  is  accomplish' 
stage  process  described  in  the  text 


The  probe  electronics  are  fairly  insensitive  to  small  changes  in  water 
conductivity  (a  10  percent  change  in  conductivity  produces  a  1  percent 
change  in  capacitance).  Furthermore,  the  apparent  absence  of  crosstalk 
between  probe  systems  makes  it  possible  to  make  precise  t M)-probe 
measurements  at  relatively  small  serrations.  In  the  past,  questions 
have  been  raised  as  to  the  accuracy  with  which  a  wire  probe  system  can 
measure  small  features  of  the  water  surface  ISturra  and  Sorrell,  1973] 
It  is  apparent  that  the  meniscus  formed  by  the  water  near  the  wire 
exhibits  a  hysteresis  which  can  couple  water  height  and  velocity  in 
such  a  way  as  to  cause  an  apparent  flattening  of  crests  and  troughs. 
This  effect  is,  however,  limited  in  magnitude  to  roughly  twice  the 
wire  diameter  and  was  therefore  discounted  as  insignificant.  In  addi¬ 
tion  it  might  be  expected  that  the  extent  of  the  meniscus  would  dis¬ 
tort  two-probe  measurements  of  wave  slope.  In  order  to  contirm  that 
these  effects  were  not  important  in  the  measurements  of  slope  pre¬ 
sented  here,  trial  measurements  were  made  at  a  variety  of  probe  spa¬ 
cing  both  greater  and  smaller  than  the  normal  1  mm.  These  measurements 
resulted  in  essentially  identical  slope  histograms,  confirming  that 
the  probe  spacing  was  substantially  greater  than  the  meniscus  effect. 

1'robe  calibration  was  accompl  ished  by  a  DC  calibration  procedure.  The 
probe  immersion  depth  in  still  water  was  varied  by  a  known  amount  and 
the  resulting  change  in  detector  output  was  gain  adjusted  to  an  ap¬ 
propriate  value.  Tills  procedure,  while  sufficient  for  simple  wave- 
height  measurements  was  found  to  be  inadequate  for  the  more  demanding 


•  differential  measurements  Involved  in  slope  analysis.  In  order  to 


calibrate  two  probe  systems  with  sufficient  accuracy  to  permit  slope 
measurement,  a  method  was  required  which  would  result  in  a  relative 
accuracy  substantially  greater  than  the  absolute  accuracy  achievable 
with  the  immersion  calibration.  The  method  used  relies  on  the  proper¬ 
ties  of  slope  statistics  in  the  following  way. 

Suppose  the  conversion  gain  of  probe  1  is  a.,  1  -  1,2,  and  the  actual 
waveheight  at  each  probe  location  is  : ^ .  The  third  moment  oi  the 
voltage  difference  at  the  detector  output  is 

hj  •  <  la^  Zx  -  a:  Z,)J  >  UM 


which  is  proportional  to  surface  slope.  Expanding  terms  one  obtains 


If  the  probe  positions  are  now  reversed,  the  same  inst rumentat ion 
gives 
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3  3  3  3  2  2 

al  <  V  >  -  a2  <  21  >  "  3al  a2  <  21  V 


+  3  aA  a2  <  Z1  Z2  > 
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which  is  equal  to  -\l  -j  only  if  aj»  a-,.  The  calibration  method  used  was 
to  calibrate  two  probes  absolutely  using  immersion  variation,  then 
fine  tune  the  gain  of  one  probe  until  the  sum  of  the  two  third  moment 
results  went  to  zero. 

Water  height  and  slope  signals  were  sampled  and  digitized  at  10  ms. 
intervals  to  12  bit  precision.  From  51,200  sample  records,  histograms 
of  surface  slope  were  constructed  to  a  resolution  of  .07  in  slope. 
These  are  given  as  Figures  24  -  33.  These  differ  substantially  from 
the  preliminary  results  included  in  the  Final  Report  of  1978  primarily 
because  of  improvements  in  instrumentation  and  redefinition  of  coor¬ 
dinates. 

Table  1  is  a  summary  of  the  mean  square  value  and  coefficient  of 
skewness  of  slope  at  a  variety  of  wind  speeds  including  those  which 
are  the  facility  standards  (those  with  "Wind  numbers").  Negative  slope 
is  defined,  rather  arbitrarily,  to  correspond  to  a  surface  rising  with 
the  wind  direction.  We  observe  that  all  such  distributions  are  nega¬ 
tively  skewed.  The  skewness  rises  with  wind  speed  up  to  a  speed  of  5 
m./sec.  and  then  falls  until  at  wind  speeds  greater  than  11.  m./sec. 
no  appreciable  skewness  is  observed.  Under  these  conditions,  the  most 
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Figure  25.  -  Histogram  of  water  wave  slopes 
taken  crosswind  at  wind  speed  5.*<  m./sec. 
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Figure  27.  -  Histogram  of  water  wave  slopes 
taken  crosswind  at  wind  speed  8.3  m./sec. 
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Figure  32.  -  Histogram  of  water  wave  slopes 
taken  downwind  at  wind  speed  8.3  m./sec. 
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Figure  33.  -  Histogram  of  water 
taken  downwind  at  wind  speed  11 


wave  slopes 
2  m./sec. 


TABLE  1:  SLOPE  DISTRIBUTION  MOMENTS 


Orientation 

Wind  No . 

Wind  Speed 

Mean  Sq . 

Skewne 

Up/Down 

I 

3.  75  m/s 

.009 

-0.4 

II 

2 

5.4 

.020 

-0.5 

II 

6.2 

.024 

II 

3 

6.7 

.028 

-0.4 

II 

6.  8 

.030 

II 

7.2 

.034 

-0.3 

II 

4 

8.3 

.040 

1 

o 

N> 

II 

9.4 

.048 

-0.2 

II 

10.7 

.057 

-0.2 

II 

5 

11.2 

.067 

-0. 1 

Cross 

1 

3.75 

.002 

II 

2 

5.4 

.004 

II 

3 

6.7 

.006 

II 

4 

CO 

co 

.012 

II 

5 

11.2 

.022 
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probable  slope  Is  slightly  positive  and  of  small  magnitude,  whereas 
large  slopes  are  more  likely  to  be  negative.  This  behavior  is 
consistent  with  results  obtained  by  Wu  [1971]  and  Bobb  et  al.  [1979] 
in  independent  model  tank  measurements.  The  largest  surface  tilt  angle 
observed  was  60°.  Mean  square  slope  increased  with  wind  speed  in  an 
almost  linear  fashion  up  to  11.  m./sec.  (see  Fig.  34).  This  is  in 
general  agreement  with  Cox  [1958]  and  Wu  [1971,1977]  although  there 
are  differences  in  absolute  magnitude  and  the  apparent  rate  of  in¬ 
crease  with  wind  speed.  The  reason  for  the  disagreement  is  unclear, 
but  could  certainly  be  due  to  differences  in  tank  construction  and 
effective  fetch.  Crosswind  slope  distributions  (Figure  34)  also  ex¬ 
hibited  an  increase  in  mean  square  slope  with  wind  speed,  but  to  much 
smaller  values. 

Sample  time  series  of  wave  height  and  slope  at  a  wind  speed  of  8.3 
m./sec.  are  given  in  Fig.  35.  The  records  are  synchronous  in  time  and 
therefore  permit  direct  comparison.  Figure  36  displays  a  similar 
comparison  with  time  derivative  of  waveheight.  These  examples  illus¬ 
trate  the  high  frequency  structure  associated  with  the  downwind  side 
of  wave  crests  which  has  been  pointed  out  by  Wu  [1977]  and  others. 
Figure  37  shows  a  synchronous  comparison  of  wave  slope  and  negative 
time  derivative.  The  similarity  of  these  records  confirms  the  notion 
that  even  small  features  of  the  surface  change  only  slowly  with  re¬ 
spect  to  their  translation  velocity  (relatively  low  dispersion)  in  the 


model  tank  situation 


Figure  35.  -  Sync 
slope  and  height 
in  the  downwind  d 
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Figure  37.  -  Synchronous  measurements  of  slope 
and  height  rate  (displayed  as  negative)  showing 
the  relative  lack  of  dispersion  at  this  scale. 
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The  statistics  described  above,  while  they  are  important  in  the  for¬ 
mulation  of  the  facet  model,  do  not  satisfy  the  requirements  of  a 
number  of  current  theoretical  formulations  predicting  Doppler  spread 
in  the  surface  scatter  channel.  These  models  generally  need  a  de¬ 
scription  of  the  surface  in  terms  of  its  second  order  space-time 
statistics.  These  are  in  some  cases  described  by  power  spectra  [Harper 
and  LaBianca,  1975]  and  in  others  by  the  space-time  correlation  func¬ 
tion  of  the  surface  [Scharf  and  Swarts,  1974].  These  statistics  con¬ 
tain  information  about  the  regularity  and  directionality  of  the  sur¬ 
face  waves  which  are  not  available  from  single  point  statistics  and 
are  therefore  necessary  in  order  to  account  for  scattering  processes 
in  which  partial  reinforcement  (Bragg)  mechanisms  are  important. 

During  the  past  year  an  extensive  set  of  measurements  was  made  of  the 
space-time  covariance  function  of  one  of  the  facility  standard  sur¬ 
faces  (Wind  4).  These  measurements  involved  the  design  and  construc¬ 
tion  of  a  computer  controlled  translator  and  software  which  permitted 
the  automatic  acquisition  of  covariances  for  128  temporal  and  32 
spatial  lags  in  one  run.  Because  of  the  relatively  low  bandwidth  of 
the  surface  motion  and  the  requirement  for  good  statistical  suffic¬ 
iency  in  the  data,  such  a  run  required  a  very  long  time  (ca.  8  hours). 
During  such  a  run,  no  other  work  could  be  done  in  the  facility. 
Therefore  the  measurements  were  taken  over  an  extended  period  of  time 
and  are  at  present  limited  to  an  ensemble  of  40,000  covariance  meas¬ 
urements  describing  a  32  cm.  radius,  90°  sector  for  time  lags  up  to 
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1.28  seconds.  This  set  of  data  describes  just  one  wind  driven  surface 
condition  and  must  be  extended  to  include  at  least  one  or  two  more. 
When  these  measurements  are  complete  they  will  be  reported  in  a  sepa¬ 
rate  technical  report. 

5.  Ultrasonic  Array  Construction:  It  has  become  apparent  that  the 
transducers  presently  available  in  the  facility  will  not  be  useful  in 
the  characterization  of  scattering  at  very  low  grazing  angles.  In 
order  to  deal  in  a  predictable  way  with  geometries  involving  low 
grazing  angles,  projector  and  receiver  beamwidths  must  be  at  least  of 
the  same  order  as  and  preferably  smaller  than  the  proposed  grazing 
angle.  During  the  past  year  the  authors  undertook  to  develop  an  ul¬ 
trasonic  array  to  serve  the  needs  of  low  grazing  angle  work  in  future 
experimentation.  This  project  was  undertaken  in  part  because  of  the 
recent  visit  to  the  facility  of  Dr.  Anthony  Pratt  of  the  University  of 
Technology,  Loughborough,  UK,  who  made  a  great  number  of  very  useful 
suggestions  based  on  his  own  experience  which  made  it  possible  for  the 
authors  to  construct  the  array  with  a  minimum  of  difficulty. 

The  array  was  constructed  using  rectangular  PZT  (5400)  transducer 
elenents  of  dimensions  23  mm.  x  9.5  mm.  and  thickness  resonant  at  320 
kHz.  The  elements  were  designed  to  be  operated  as  quarter  wave  air- 
backed  resonators.  The  length  and  width  were  chosen  so  as  suppress  any 
other  resonance  modes  at  operating  frequency.  The  transducers  are 
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mounted  in  a  foam  filled  polycarbonate  (Lexan)  housing.  Although  the 
acoustic  impedance  of  the  foam  backing  is  much  higher  than  that  of 
air,  it  is  still  an  order  of  magnitude  smaller  than  that  of  water  and 
thus  suffices  as  a  backing  material. 

The  construction  details  of  the  array  are  shown  in  Figure  38.  Hie 
housing  is  built  of  3/8"  thick  Lexan  plate.  It  is  a  partitioned  rec¬ 
tangular  box  with  both  sides,  open  filled  with  a  self-foaming  resin, 
milled  out  to  accommodate  the  transducer  elements.  The  elements  were 
mounted  on  1  cm.  centers  in  the  face  of  the  housing  and  the  front  face 
was  sealed  with  a  thin  layer  of  silicone  rubber  (RTV).  Electrical 
connection  was  made  by  soldered  connections  carried  to  a  cable  assem¬ 
bly  of  17  -  RG174A  coaxial  cables.  The  array  is  amplitude  shaded  by  an 
assembly  of  seventeen  precision  potentiometers  (10  turn  Helipot) 
connected  directly  to  the  cable  assembly. 

The  array  was  characterized  using  sinusoidal  bursts  about  2  ms.  long 
at  100  ms.  intervals.  All  measurements  were  taken  with  the  receiver 
about  3  m.  distant  from  the  array.  After  individual  impedance  meas¬ 
urements  were  made  to  verify  the  manageability  of  the  design,  indi¬ 
vidual  element  frequency  responses  were  taken  in  the  finished  array.  A 
selection  of  these  is  given  as  Figure  39.  Using  this  information  the 
shader  was  adjusted  to  give  an  effectively  uniform  shading  function  at 
the  resonant  frequency.  The  frequency  response  of  the  array  on  axis 
under  these  conditions  is  given  as  Figure  40. 
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Figure  38.  -  Details  of  construction  of  the  17 
element  320  kHz.  array.  The  transducer  elements 
are  PZT  (5400). 
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Figure  39.  -  Responses  of  a  selection  of  five 
transducer  elements.  Some  portion  of  the  non¬ 
uniformity  could  have  been  due  to  construction 
imperfections. 
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Figure  40.  -  Overall  response  of  the  finished 
array  with  all  elements  driven  so  as  to  generate 
equal  output.. 
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Figure  41.  gives  the  polar  response  at  resonance  along  the  axis  of  the 
array  of  a  single  element.  The  -10  dB.  beam  width  is  seen  to  be  35°. 
Figure  42.  gives  the  polar  response  of  the  uniformly  shaded  array  at 
322.6  kHz.  Under  these  conditions  the  main  lobe  -10  dB.  beam  width  is 
5°,  but  there  are  sidelobes  at  29°  at  only  5  dB.  below  the  main  lobe. 
In  addition  the  power  radiated  between  the  principal  lobes  is  only  17 
dB.  down.  The  pattern  perpendicular  to  the  principal  axis  of  the  of 
the  array  is  shown  in  Figure  43. 

Several  experiments  with  shading  were  made  to  test  the  effect  of 
shading  on  the  pattern  of  the  array.  The  first  of  these  was  a  cosine 
shading  function  used  as  a  basis  for  comparison  with  more  complicated 
shading  functions.  The  result  is  given  in  Figure  44.  As  can  be  seen 
the  simple  shading  function  resulted  in  a  reduction  of  main  lobe  width 
to  4°  and  a  substantial  reduction  in  the  near  sidelobe  power  to  at 
most  -  21  dB.  A  second  attempt  at  shading  was  the  application  of  the 
Dolf-Chebychev  method  for  a  predicted  sidelobe  rejection  of  30  dB.  The 
result  of  this  experiment  is  given  as  Figure  45.  The  actual  reduction 
obtained  was  only  -23  dB.,  possibly  because  of  approximations  made  in 
adjusting  the  shading  potentiometers.  Again  the  mainlobe  width  was 
about  4°. 

The  general  result  of  this  adventure  in  array  design  was  the  acquisi¬ 
tion  of  a  new  and  useful  tool  for  narrow  beam  probing  of  scattering 
surfaces.  The  array  appears  to  have  sufficient  efficiency  and  band- 
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width  for  pulsed  monochromatic  probing  of  the  surface  at  low  grazing 
angles.  Not  surprisingly,  the  patterns  achieved  with  shading  do  not 
appear  terribly  sensitive  to  the  choice  of  shading  function. 

6.  Frequency  Spreading  Studies:  Rather  independently  of  the  above 
work,  previous  studies  of  frequency  spreading  in  forward  scatter  have 
been  continued  in  order  to  explore  an  anomalous  shift  in  mean  received 
frequency  under  conditions  of  geometric  symmetry  [Zornig,  1978].  These 
results  are  as  yet  not  fully  explained.  However,  a  theory  explaining 
some  of  the  effects  observed  is  in  the  process  of  publication  [Tuteur 
et  al . ,  1979]  and  will  be  available  in  the  near  future. 

The  experiment  conducted  was  basically  the  same  as  that  described  in 
the  above  referenced  work.  A  preemphasized  probing  signal  s(t)  is 
designed  and  synthesized.  This  signal  is  applied  to  a  projector, 
interacts  with  the  surface,  is  received  by  a  hydrophone  and  is  digi¬ 
tized  at  the  sampling  rate  used  in  its  synthesis.  Because  of  the 
linearity  of  the  process,  it  is  assumed  and  has  been  demonstrated 
[Zornig  and  McDonald,  1974]  that  the  signal  waveform  received  after 
digitization  is  equal  to  the  convolution  of  the  original  probing 
signal,  the  impulse  responses  of  the  two  transducers,  and  the  time 
varying  impulse  response  of  the  surface.  Since  the  surface  impulse 
response  is  slowly  varying  with  respect  to  the  time  scale  of  the 
signal  itself,  it  is  possible  to  represent  this  process  after  trans- 
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formation  to  the  signal  frequency  domain  as  follows: 
R(un,t)  -  s(«n)  Pq(“)  H (u> , t )  P1(ui) 


2H 
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where  S(wn)  is  the  probing  signal  applied  to  the  projector 
Pq(w)  is  the  projector  transfer  function 
Pj(w)  is  the  receiver  transfer  function 

H(u,t)  is  the  time  varying  transfer  function  of  the  surface 
RC^nif)  Is  the  signal  out  of  the  hydrophone  and  amplifiers 

The  data  collection  instrumentation  removed  a  fixed  time  delay  from 
the  received  signals  which  is  reflected  in  the  notation  here  in  the 
use  of  the  same  t  in  both  H(w,t)  and  R(wn,t)  and  the  suppression  of 
time  dependence  in  S(wn) .  The  representation  of  S(uin)  and  R(u>n,t)  as 
defined  on  a  finite  domain  of  frequencies  is  based  on  the  fact  that 
both  waveforms  are  wholly  contained  in  relatively  short  time  windows 
and  are  analyzed  in  terms  of  discrete  Fourier  transforms. 

In  this  representation,  H(w,t)  is  a  stationary  ergodlc  random  process 
expressing  the  amplitude  and  phase  response  of  the  surface  to  fre¬ 
quency  w  at  time  t.  Care  was  taken,  of  course,  to  sample  in  t  at  rates 
in  excess  of  the  Nyquist  rate  of  the  process;  the  actual  sampling  rate 
used  varied  between  50  and  300  per  second.  In  practice,  sampling  in  t 
was  limited  to  64  contiguous  samples  in  a  burst.  This  limit  was  im¬ 
posed  by  memory  availability  in  the  data  acquisition  computer,  and  is 
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equivalent  to  CW  pulse  lengths  of  from  850  to  20,000  cycles  of  the 
carrier.  The  recorded  samples  of  H(w,t)  were  then  Fourier  transformed 
on  the  variable  t  to  give  an  ensemble  of  instantaneous  spreading 
functions. 

4(VV  “  <  F(vV  *  W(V  *  (20) 

f 

Ensembles  of  roughly  200  measurements  of  were  then  windowed 

and  averaged  to  give  the  frequency  spreading  function, 

N-l  -jv.mAt 

F(«_,v  )  -  (S(u  )P_(u  )  P.(u  ))_i  •  (dt  2  R(«  ,t  )e  k  )  (19) 

n  k  n  u  n  l  n  m=o  n  m 

In  previous  work,  a  systematic  shift  of  the  frequency  of  the  maximum 
of  the  frequency  spreading  function  had  been  observed  under  some 
symmetric  geometries.  Since  a  number  of  theoretical  formulations 
predicted  that  under  conditions  of  symmetric  transducer  location 
Doppler  spread  would  be  symmetric,  it  was  suggested  that  the  observed 
asymmetries  could  be  due  to  small  misalignments  of  the  experimental 
apparatus.  Therefore  a  careful  set  of  experiments  was  performed  in 
which  transducer  alignment  was  surveyed  to  a  high  degree  of  precision 
and  then  varied  so  as  to  test  the  sensitivity  of  the  results  to 
alignment  error.  Figure  46.  gives  a  set  of  Doppler  spread  spectra 
taken  crosswind  at  4  harmonically  related  frequencies  and  relatively 
high  roughnesses.  The  apparatus  was  exactly  aligned  both  with  with 
respect  to  itself  and  to  the  wind  direction.  It  can  be  seen  that  a 
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Figure  47(a):  -  Frequency  spread  spectrum 
taken  at  17  deg.  grazing  angle  and  a 
6  deg.  "dogleg"  alignment  shift  from  a 
crosswind  orientation. 
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Figure  47(b):  -  Frequency  spread  spectrum 
taken  at  17  deg.  grazing  angle  and  a  6  deg. 
"dogleg"  alignment  shift  from  a  crosswind 
orientation . 
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Figure  *»9:  -  Mean  frequency  shift  and  RMS 
spread  at  17  deg.  grazing  angle  for  several 
counterclockwise  "dogleg"  alignment  shifts 
from  a  crosswind  orientation. 


small  up  shift  Is  apparent,  Increasing  with  carrier  frequency. 

In  order  to  test  sensitivity  to  alignment,  the  receiver  was  moved  6° 
in  azimuth  while  maintaining  grazing  angle  and  range.  This  type  of 
misalignment  has  the  effect  of  shifting  the  center  of  the  Fresnel 
zones  on  the  surface,  the  specular  point,  out  of  the  center  of  the 
transducer  beam  patterns  and  of  shifting  the  specular  path  with  re¬ 
spect  to  the  wind  direction,  but  only  slightly.  The  results,  shown  in 
Figure  47,  demonstrate  the  very  large  change  in  Doppler  shift  produced 
by  a  misalignment  of  this  type.  When  carefully  surveyed,  the  transdu¬ 
cers  are  oriented  to  within  a  small  fraction  of  a  degree.  A  summary  of 
this  experiment  repeated  for  several  "dogleg"  angles  is  given  as 
Figures  48  -  49.  As  can  be  seen  the  results  are  quite  regular. 

In  order  to  test  the  sensitivity  of  the  spreading  to  wind  direction  a 
number  of  runs  were  taken  at  specular  alignment,  but  at  various 
orientations  to  the  wind.  Several  of  these  are  given  as  Figure  50  for 
a  single  carrier  frequency.  As  can  be  seen,  the  neither  the  spread  nor 
the  shift  is  sensitive  to  the  exact  wind  direction. 

Finally,  the  transducer  alignment  experiment  was  repeated  in  the 
up/down  wind  orientation.  The  result,  Figures  51  -  55,  showed  only  a 
very  small  sensitivity. 

The  particular  sensitivity  in  the  crosswind  direction,  the  authors 
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Figure  51(b):  -  Frequency  spread  spectrum 
taken  at  17  deg.  grazing  angle  and  well 
surveyed  downwind  orientation. 
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Figure  52(b):  -  Frequency  spread  spectrum 
taken  at  17  deg.  grazing  angle  and  a  6  deg. 
"dogleg"  alignment  shift  from  a  downwind 
orientation. 


SPREAD  (HZ)  SHIFT  (HZ) 

5-  50-  75-  100-  -20-0  0-0  ! 


SPREAD  (HZ) 


0*25  0-50  0*75  1-00 

INCIDENT  FREQ-  (MHZ) 


- ! - 1 - 1 - 1 - 

0-25  0-50  0*75  1-00 

INCIDENT  FREQ-  (MHZ) 

Figure  51*:  -  Mean  frequency  shift  and  RMS 
spread  at  17  deg.  grazing  angle  for  several 
counterclockwise  "dogleg"  alignment  shifts 
from  a  downwind  orientation. 
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Figure  55:  -  Comparison  of  frequency  shift 
ana  spread  observed  at  exact  transducer 
alignment  in  crosswind  and  up_/downwind 
orientation. 
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believe,  can  only  be  attributed  to  asymmetry  in  the  wind  driven  sur¬ 
face.  The  exact  mechanism  is  the  subject  of  the  theoretical  investi¬ 
gation  referred  to  above.  In  general  qualitative  terms  the  shift  in 
the  case  of  exact  alignment  can  be  seen  as  the  result  of  unequal 
contribution  to  the  scattered  signal  from  portions  of  the  surface  up 
and  down  wind  from  the  specular  point.  This  unequal  contribution,  we 
suggest,  is  due  to  the  asymmetry  of  the  shape  of  the  surface  that 
result  from  the  fact  that  it  is  wind  stressed.  In  the  case  of  mis¬ 
alignment,  the  shift  probably  results  from  asymmetric  weighting  of 
contributions  from  the  surface  due  to  the  transducer  pattern  mis¬ 
alignment;  wind  direction  has  apparently  little  effect. 

7.  Conclusion:  It  can  readily  be  seen  that  this  report  is  a  summary  of 
the  current  state  of  a  number  of  projects  in  progress.  The  main  thrust 
of  the  contract  under  which  this  work  has  been  supported  has  been  a 
better  understanding  of  the  reverberation  from  rough  surfaces,  espe¬ 
cially  in  bistatic  configurations  and  at  low  grazing  angles.  The 
scattering  strength  measurements,  surface  slope  measurements  and  facet 
model  development  have  led  us  to  believe  that  the  primary  mechanism 
involved  in  low  angle  backscatter  is  diffraction  and  that  the  model 
has  a  high  probability  of  developing  into  a  useful  predictive  tool. 
However,  at  present  there  are  substantial  analytical  obstacles  to  the 
development  of  a  full  facet  model  which  must  either  be  overcome  or 
sidestepped.  It  seems  clear  that  the  scattering  strength  observed  at 
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long  wavelengths  has  turned  out  to  be  considerably  higher  than  that 
measured  at  very  short  wavelengths.  However,  additional  measurements 
at  a  variety  of  wavelengths  both  intermediate  and  much  longer  are 
feasible  and  desirable  as  guidance  in  testing  the  model. 

The  nature  of  frequency  spreading  behavior  has  become  somewhat  clearer 
as  a  result  of  the  sensitivity  tests,  but  a  full  explanation  is  still 
not  available.  Recently,  two  ocean  experiments  have  been  conducted 
which  will  provide  an  additional  basis  for  comparison  with  theory.  In 
addition,  the  availability  of  a  complete  second  order  description  of 
the  model  tank  surfaces  used  in  these  studies  will  permit  comparison 
with  general  Helmholtz  integral  formulations.  The  authors  wish  to 
point  out  that  the  data  presented  in  this  report  are  representative  of 
a  large  base  of  data  obtained  through  experimentation  and  analysis 
which  has  not  been  included  for  reasons  of  clarity  and  economy.  These 
data  are,  however,  available  to  interested  researchers  on  request. 
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